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2.1 Introduction 



Thermodynamics is an efficient means of finding most general regularities in the theory of 
thermo electricity. 

In 1854, W. Thomson using equilibrium thermodynamics established the interrelation between three 
thermoelectric effects that arise in a thermocouple thermoelectric circuit with the legs of individual 
isotropic media (Figure 2.1) in the form of relationships 

dWdT - a- t 
dn/dT - II/T-1 











(2.1) 
(2.2) 



where a is the Seebeck coefficient, 77 is the Peltier coefficient, and Tis the Thomson coefficient. Equation 
2.1 and Equation 2.2 imply two relationships of frequent use in thermoelectricity: 



77 = aT is the first Thomson relationship 
t = T(da/dT) is the second Thomson relationship 



(2.3) 
(2.4) 



It is common knowledge today that a rigorous thermodynamic theory of thermoelectricity can be 
constructed only on the basis of a nonclassical, more general theory of macroscopic description of 
nonequilibrium processes called nonequilibrium thermodynamics or thermodynamics of irreversible 
processes. 2-7 

Stated below are the basic concepts of modern nonequilibrium thermodynamics as applied to 
thermoelectricity, which make possible specification and generalization of the Thomson theory not only 



2-1 
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FIGURE 2.1 Thermocouple thermoelectric circuit with the legs of individual isotropic media A and B. 
Thermocouple junctions are at temperatures T and T+dT. 

for thermocouple power converters, but also for thermoelectric media of any complexity used in them, 
generalization of fundamental thermoelectric relationships, discovery of new thermoelectric effects, and 
formulation of general efficiency criteria for materials used in thermoelectric devices. 

2.2 Nonequilibrium Thermodynamics: General Concepts 

This is a reminder in brief the concepts of classical thermodynamics. 

The first law of thermodynamics is a concept expressing the law of conservation of energy whereby the 
quantity of system-absorbed heat 8Q can be expressed as the relationship 

8Q = dU + 8A - fjudn (2.5) 

where dU is a change in internal system energy, 8A is work performed by the system against external 
forces, /a is system chemical potential, and d« is a change in the number of particles in the system. 
The second law of thermodynamics, whereby 

TdS > SQ (2.6) 

where the equality sign is related to reversible, and nonequality sign, to irreversible processes; T is 
absolute temperature, dS is total differential of state function S which is commonly called entropy. The 
classical theory of thermodynamics of reversible processes outlines only the idealized processes. For real, 
i.e., nonequilibrium processes (such as heat transfer, electric current flow), classical theory provides only 
the inequalities of 2.6 type, indicating the possible evolution of these processes. 

A decisive step for the creation of nonequilibrium thermodynamics was made in 1931 by Onsager. 
He formulated for thermodynamics a new general principle of minimum energy dissipation. Besides, 
Onsager put forward a partial principle that was a generalization of Ohm's and Fourier equations. 
According to this partial principle, irreversible processes can be described by linear differential equations 
with constant coefficients, L ik , for which the symmetry relations of the type L ik = L ki are satisfied. The 
Onsager principle was of considerable importance in the development of thermodynamics of irreversible 
processes. 

Unlike Onsager, Prigozhin 4 put forward a new, also general principle for thermodynamics, namely the 
principle of minimum entropy production, which proved to be more convenient for the solution of 
practical tasks than the Onsager principle. This principle determines a criterion for thermodynamic 
system evolution whereby under nonequilibrium conditions held constant a thermodynamic system 
evolves to a steady-state having relative stability and characterized by a minimum rate of entropy 
production. The most essential conclusion in the Prigozhin theory is that with increasing the degree of 
system state nonequilibrium, for example, with increasing thermodynamic forces, the state of 
homogeneous chaos becomes unsteady, kinetic-phase transitions occur in the system, and ordered 
dissipation structures are formed providing a higher entropy production rate as compared to the chaotic 
state. Numerous investigations show that this regularity is universal, valid for various-complexity 

9,10 

systems, 
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At the present time linear and nonlinear versions are recognized in nonequilibrium thermodynamics. 
The best developed is a linear nonequilibrium thermodynamics with its prime objective to study 
nonequilibrium processes for near-equilibrium states. 

In thermodynamics of irreversible processes the systems with nonequilibrium processes are regarded 
as continuous media, and their state parameters as field variables, i.e., continuous coordinate and time 
functions. 

For a macroscopic description of nonequilibrium processes the following method is used: a system is 
represented as consisting of elementary volumes (Figure 2.2), although big enough to comprise a very 
large number of particles. The state of each elementary volume in a medium is characterized by 
temperature, pressure, and other thermodynamic parameters depending on coordinates and time, i.e., it 
is assumed that this element is actually in the state of local equilibrium even when the system in general is 
thermally nonequilibrium. This assumption is called the principle of local equilibrium. 

Quantitatively, the nonequilibrium processes in this method are described by balance equations for 
elementary volumes based on the laws of conservation of mass, momentum, and energy, entropy balance 
equations, as well as phenomenological equations of the processes in question that express mass, 
momentum, electric current, and energy flows in terms of the gradients of thermodynamic parameters. 

The methods of theory of irreversible processes allow formulating for nonequilibrium processes the first 
and second laws of thermodynamics in a local form, as well as obtaining from the general principles a 
complete set of energy, mass, momentum, and entropy transport equations, i.e., equations for the description 
of thermoelectric phenomena, hydrodynamics, thermal conductivity, diffusion, etc. These equations are 
fundamental for nonequilibrium thermodynamics. In thermoelectricity, where mass and momentum 
transport is not essential, local equations of energy balance and entropy are of basic importance. 

The law of conservation of energy for volume elements is the first law of thermodynamics in the theory of 
irreversible processes. Here one should take into account that total unit energy is composed of unit 
kinetic, unit potential energy in the field of forces F*, unit internal energy u, which is the energy of 
thermal motion of particles, and the average interaction energy. Between them the balance equation for u 
is of the form 

p(dpu/df) = -div(p«v + q) - £ P afs (dv a /dxp) + X JA (2-7) 

a/3 k 

where q is density, v is mass transport velocity, and in so doing, the rate of change in momentum density 
per one particle dpu/dt is determined by the divergence of internal energy flows pux and heat flow q, and 
the work of internal mechanical stresses Y.a/3 ^as(dv a /9Xg) and external forces Xk JA- 



T+dT 




FIGURE 2.2 Elementary volume dV of thermoelectric medium considered in nonequilibrium thermodynamics. 
This volume is also used in local thermoelement models. 
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Entropy balance equation. In thermodynamics of irreversible processes it is assumed that the entropy of 
volume unit s (local entropy) is the same function of internal energy u, unit volume v — lip and 
concentration c, as in the state of full equilibrium; hence, regular equalities of classical thermodynamics 
hold true for it. These concepts together with the laws of conservation of mass, momentum, and energy 
permit finding the entropy balance equation 

p(ds s /dt) = -div J s + a s (2.8) 

where a s is unit local entropy production per unit time, J s is entropy flow density expressed in terms 
of heat flow density, diffusion density, and that part of stress tensor which is related to nonequilibrium 
processes (i.e., in terms of elastic stress tensor 77 a/3 ). Entropy (in contrast to mass, energy, and 
momentum) is not retained, but increases with time in volume element due to irreversible processes with 
the rate of cr s increase. Besides, entropy can change on its leaking into or out of volume element, which is 
not related to irreversible processes. The positive character of entropy production cr < is expressed in 
thermodynamics of irreversible processes by the law of entropy increase. 

Entropy production cr s is determined only by irreversible processes (such as diffusion, thermal 
conductivity, viscosity) and is equal to 



iJ.x, 



(2.9) 



where J, are flows (for example, diffusion flow \, heat flow q, elastic stresses tensor TI a p), and X ; are their 
conjugate thermodynamic forces, i.e., gradients of thermodynamic parameters causing a deviation from 
equilibrium state. To obtain in thermodynamics of irreversible processes a closed set of equations 
describing nonequilibrium processes, the flows of physical quantities should be expressed by 
phenomenological equations in terms of thermodynamic forces. 

Phenomenological equations. Thermodynamics of irreversible processes proceeds on the assumption 
that with low system deviations from thermodynamic equilibrium the arising flows are linearly 
dependent on thermodynamic forces and described by phenomenological equations of the type 

J, = XL,A (2.10) 

k 

where L ik are kinetic (phenomenological) or transport coefficients. In direct processes thermodynamic 
force X^ brings about flow J^., for example, temperature gradient brings about heat flow (thermal 
conductivity), concentration gradient — substance flow (diffusion), velocity gradient — momentum 
flow (that determines viscosity), electric field — electric current (electric conductivity). Such processes 
are characterized by kinetic coefficients, proportional to coefficients of thermal conductivity, diffusion, 
viscosity, and electric conductivity. These coefficients are also called kinetic coefficients or transport 
coefficients. Thermodynamic force X k can induce thermoelectric current, and electrochemical potential 
gradient — heat flow (Peltier effect). Such processes are called cross-flow or overlapping effects. They are 
characterized by coefficient L ik with i ¥^ k. With regard for phenomenological equations, entropy 
production is equal to 

CTs = 2>,I,A > (2.11) 

i,k 

In steady-state the value cr s is minimum under assigned external conditions preventing from establishment 
of equilibrium. This statement is known as Prigozhin's theorem. In the state of thermodynamic 
equilibrium, cr s — 0. In the above- discussed examples thermodynamic parameters are continuous 
coordinate functions. Nonequilibrium systems are possible where thermodynamic parameters are 
changed abruptly (heterogeneous systems), for example, electron gas in metal thermocouple or gases in 
vessels connected with a capillary or membrane. If temperatures T and chemical potentials jjl of 
components on both sides of the thermocouple junction or membrane between the vessels are not equal 
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(T l > T 2 and /jl 1 > jj, 2 ), thermodynamic forces (X n = \IT 2 ~ l/T 1 , X m = fi 2 /T 2 — /jl 1 /T 1 ) bring about 
particle and energy flows (J m = L 11 X,„+L 12 X„, J„ =L 21 K m +L 22 K n ) between the heterogeneous system 
parts, creating thermal difference in pressures. In these examples the flows and thermodynamic forces are 
scalars and accordingly such processes are called scalar. In thermoelectricity the flows and thermodynamic 
forces are vectors; hence they are called vector processes. In hydrodynamics at shear viscosity 
thermodynamic forces and flows are tensors. 

Thermodynamics of irreversible processes provides a theoretical basis for the investigation of 
open systems. It offers explanation of many nonequilibrium phenomena in conductors, including 
thermoelectric, galvanomagnetic, and galvanothermomagnetic phenomena. 

According to the thermodynamics of irreversible processes, 2-6 the basic equation of thermodynamics 
for quasistatic processes is applicable to parts of physical system (medium) being in local equilibrium 

TdS = dU + 8A - £ i±,.dn r (2.12) 

r 

that follows from Equation 2.5 and Equation 2.6. 

Thermodynamic flow caused by thermodynamic forces X k is related to them by a linear law 
(Equation 2.10) where kinetic coefficients L ik describing medium properties are assumed to be known. 

It should be noted that even in modern publications, a linear law of type shown in Equation 2.10 is 
not infrequently written in a simplified form, suitable only for isotropic medium. However, for a full and 
correct description it should be assumed that thermodynamic flows and forces are components of three- 
dimensional vectors, and kinetic coefficients — components of second-rank tensors so that summation 
limit in Equation 2.10 is equal to the number of all vector components of thermodynamic forces. 

If thermodynamic flows and forces are selected such that entropy production rate cr s can be 
represented as Equation 2.11, for kinetic coefficients the Onsager symmetry principle holds true: 

L ik = L ki (2.13) 

In the presence of a magnetic field bold dependence of kinetic coefficients more general relations hold 
true: 

I*(B) = L fa -(-B) (2.14) 

where B is magnetic induction vector. 

As a rule, a relation between electric field E, electric current density j, and conducting medium 
property (electric resistivity p), as well as a relation between temperature gradient VT, heat flow density q, 
and medium thermal conductivity k is found from Ohm's and Fourier laws, which are of the form 

E = pj (2.15) 

q=-KVT (2.16) 

Equation 2.15 and Equation 2.16 are valid for the cases when either electric processes or thermal 
processes alone occur in a medium; therefore these equations in principle cannot describe thermoelectric 
effects for which the simultaneous existence and interaction of thermal and electric processes is 
important. 

2.3 Thermodynamic Theory of Thermoelectricity 
2.3.1 Generalized Ohm's and Fourier Laws 

Consider generalizations of Equation 2.15 and Equation 2.16 for thermoelectric phenomena. 

The application of formula 2.15 and formula 2.16 to heat and electric current dissipation processes 
yields the following basic equations of thermoelectricity 12-14 : 

Ej = Pikik + a im (dT/dx m ), q t = II lk j k - fc /m (8T/9x m ) (2.17) 
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The first of relations given in Equation 2.17 is often written as 

h = ViA - o- ik a km (dT/dx„) (2.18) 

In Equation 2.17 and Equation 2.18, E; is component of electric field intensity, p ik are components of 
electric resistivity tensor, Oft are components of electric conductivity tensor, qi is the component of heat 
flow density vector, j k is the component of electric current density vector, K im is the component of 
thermal conductivity tensor, x m are Cartesian coordinates, indexes i,k,l,m run the values 1, 2, 3, the 
summation being done with respect to doubly encountered indexes. 

Relations given in Equation 2.17 are called generalized Ohm's and Fourier laws. 

As compared to conventional Ohm's law (Equation 2.15) and Fourier law (Equation 2.16), the 
generalized laws (Equation 2.17) include additional components comprising coefficients a ik and H ik . 
These are the Seebeck and Peltier coefficients. Their physical meaning and relation to other kinetic 
coefficients will be considered below. 

The Onsager principle (Equation 2.13) implies the following properties of kinetic coefficient tensors p, 
k, a, a, 77: 

ojt(B) = ofe(-B), p ik (B) = Pfc(-B), Klm (B) = k,„,(-B), a km (B) = (VT)II mk (-B) (2.19) 

For B = formulae 2.19 yield the following symmetry relations for tensors p and /?: 

Pik = Pkh «lm = "ml (2-20) 

However, it should be borne in mind that tensor a, generally speaking, is not symmetrical even in the 
absence of a magnetic field. 

2.3.2 Generalization of the Thomson Relations 

The Thomson relations establish fundamental laws of relation between the Seebeck, Peltier, and 
Thomson effects. These relations are widely used in the theory of thermoelements and in the calculation 
of their parameters (efficiency, coefficient of performance, irreversible losses through caloric 
thermoelectric effects, etc.). 

The formulation of the Thomson relations has experienced essential evolution. These relations were 
first established by Thomson based on the use of classical thermodynamics of equilibrium reversible 
processes and looked at as a relation between the Peltier and Seebeck coefficients (the first Thomson 
relation [Equation 2.3]) and a relation between the Thomson and Seebeck coefficients (the second 
Thomson relation [Equation 2.4]) for isotropic media. 

As indicated above, classical thermodynamics of reversible equilibrium processes is not quite correct 
for thermoelectric phenomena, so it was natural to study the Thomson relations based on the 
thermodynamics of irreversible processes after its fundamentals had been created by Onsager. This work 
was performed by Domenicali and Samoilovich and Korenblit [12]. Nevertheless, the employment of 
the theory of irreversible processes did not change the expression of the Thomson relations (Equation 2.3 
and Equation 2.4). 

The reason for this is that the Seebeck, Peltier, and Thomson effects as such, although always 
accompanied by nonequilibrium electric conductivity and thermal conductivity processes, in the 
isotropic medium in a steady-state can be regarded as "quasi-equilibrium" and "quasireversible" in the 
sense that under these conditions their contribution to local entropy change s of unit volume is zero. 

On substituting into Equation 2.21 heat flow vector from Equation 2.17 
we get 

s = -div kV In T + pf/T (2.21) 

From Equation 2.22 it can be seen that in this case the unit change of entropy is determined only 
by two nonthermoelectric medium parameters, i.e., thermal conductivity and electric resistivity. 
Therefore, if we consider by convention only thermoelectric phenomena, abstracting from the 
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accompanying nonequilibrium processes of electric conductivity and thermal conductivity, the second 
law of thermodynamics for thermoelectric components of processes under study can be used as a strict 
equality ds = — divq/T, valid for equilibrium processes instead of inequality d5 > — div q/T of 
nonequilibrium thermodynamics. 

Thus, a deeper insight into the Thomson relations in the first stage of thermodynamics of irreversible 
processes application to thermoelectricity did not culminate in the specification or establishment of new 
regularities. 

The authors of Ref. [13] formulated the Thomson relations for the anisotropic media and the presence 
of a magnetic field. The relations generalized for this case took the form 

77(B) = rS(-B) (2.22) 

f r (B) = T[M(-B)/bT] (2.23) 

where the tilde sign signifies tensor transposition. 

Another generalization of the Thomson relations, different from Equation 2.22 and Equation 2.23, was 
proposed in Ref. [14]. The necessity and content of this generalization lies in the following. 

As can be seen from Equation 2.22 and Equation 2.23, the right-hand and left-hand tensors in these 
relations differ in the transposition and sign of vector B. Therefore, physically Equation 2.22 and 
Equation 2.23 are formal statements on the equality of values belonging to two absolutely different 
physical states of medium. The first state (left-hand side of equalities) is characterized by a certain 
orientation of external influence vectors j, q (or VT) and B, superimposed on the medium. The second 
state (right-hand side of equalities) is different in that vector j here is oriented as vector q of the first state; 
vector q here is oriented as vector j of the first state, and magnetic field B reversed its direction. Thus, 
directly from Equation 2.22 and Equation 2.23, a relationship between the Seebeck, Peltier, and Thomson 
coefficients for one certain physical state of medium is not established, unlike the original Thomson 
formula given in Equation 2.3 and Equation 2.4. 

Therefore, let us introduce into our consideration the local Peltier, Seebeck, Thomson, electric 
resistivity, and thermal conductivity coefficients for an arbitrary, but quite certain orientation of vectors j, 
VT, q, and B and reformulate the Thomson relations in these terms. These coefficients are of the form 

P = Jo(P jo), k* = t (£ t ), t = jo(Tt ), a* = t (a j ), 77* = t (77j ) (2.24) 

where t and j are unit vectors of temperature gradient and electric current density. 
Let us find a new form of the first Thomson relation. 
Let us consider coefficient 77* to the product a "T with regard for Equation 2.22: 

a*T= rt (aj ) = 77*S (2.25) 

where 

8= a*(B)/a*(-B) (2.26) 

A dimensionless parameter 5 can be represented as 

S(j , to, B) = j a(B)t /j a(-B)t (2.27) 

whence it is seen that this parameter in the absence of a magnetic field 

8(j ,T ,0)=l (2.28) 

In the presence of a magnetic field 5 can be a value different from unity only in the case when 

a(B) 7^ a(-B) (2.29) 
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This is possible if a magnetic field reversal causes a change in the Seebeck coefficients (Umkehr effect in 
low-symmetry crystals), whereby 

-1<5<1 (2.30) 

from (3.10). With regard for thermodynamic limits for S (3.15) it is evident that 

n* > a*T (2.31) 

Thus, in the general case in the presence of a magnetic field in low-symmetry crystals a new form of 
relation (Equation 2.31), in contrast to Equation 2.3, takes the form of inequality. Only for the isotropic 
medium Equation 2.31 takes the form of absolute Thomson equality (Equation 2.3). 

Equation 2.31 is a new form of the first Thomson relation. 

Let us use the above approach for the second Thomson relation as well. 

Dividing the Thomson and Seebeck tensors into symmetrical and nonsymmetrical, as well as even and 
odd components, from Equation 2.23 we get 

t* = T(da* + /dT) - T(da~/dT) - 2a*~ (2.32) 

where 

a + = t (d + j ), a~ = to(cT j ) (2.33) 

Equation 2.33 is a new form of the second Thomson relation. 

In contrast to Equation 2.23, it comprises two odd with respect to magnetic field components 
— T(da"~/dT) — 2a*~ . In so doing, the presence of — 2a*~ component allows for the existence of an 
odd Thomson effect in a medium where thermoemf is temperature independent. 

Thus, a new form of the Thomson relations yields two essential results, namely the inequality 
77* > a*T and the Thomson effect odd with respect to a magnetic field. 

2.3.3 Thermodynamics of Caloric Effects in Thermoelectricity 

The value of heat release rate Q in the anisotropic medium can be represented as 

Qj = -div(w) - div kVT (2.34) 

where energy flow vector is determined by the expression 

w = - kVT + ft j + (/Z/e)j (2.35) 

/I is electrochemical potential of current carriers, e is their charge, and k is the thermal conductivity 
tensor. Expression 2.34 with regard for Equation 2.35 can be rearranged as 

Qj = On + Qt + Qb + Qy (2.36) 

where Q[j is the Peltier heat, Q T is the Thomson heat, Q B is the Bridgman heat, and Qj is the Joule heat 
determined by the relations 

Qn = -rj{Div r (a + - a 5- ) - rot(N + - IsT)} (2.37) 

Q T = -Tj{rXa s+ - cO/3T-2a s ~}VT- T{3(N + - N~)3T - 2N~}[Vrx j] (2.38) 

Q B = -T Dev(a + - d s ")/Def j - T(N + - N~)rot j (2.39) 

Qy = )P S ) (2-40) 

i 2006 by Taylor & Francis Group, LLC 



Modern Thermodynamic Theory of Thermoelectricity 



2-9 



TABLE 2.1 Caloric Thermoelectric Effects 



Name 



Effect Type 



Specific Heat Production 



Peltier effect 






Divergent 
Eddy 






Even 
Odd 
Even 
Odd 






-TjDiv T a s+ 
TjDiv T a'~ 
TjrotN + 
-TjrotN" 


Bridgman 






Deviation 
Eddy 

Longitudinal 






Even 
Odd 
Even 
Odd 

Even 






-TDev(a s+ ) :Defj 
TDev(a s ~) : Def j 
-TN+rotj 
TTSTrotj 

'dT 


Thomson effect 












Odd 




















Even 












Transverse 
























Odd 






./3« s ., \ 


In Equation 


2.36 


through Equation 


2.40, 


(D 


iv r t)-, 


= 01 


ij/dXj) 


r=const. (Devf s ), ; - 



dT 



N + [r-j] 



») 



; - (V3)s P f s , 

(Def r ?),%• = l/2(9/;/3x- + djJdXi), colon signifies a biscalar product of tensors. The formulae for heat 
release rate are summarized in Table 2.1. From the table it is evident that essential importance for caloric 
effects is assumed by spatial derivatives of these vectors, as well as by gradients of tensor components a 5 ~ , 
resulting in new effects missing from classical theory. 



2.3.4 Local and Integral Efficiency of Thermoelements 

There are two approaches with two models developed for the description of thermoelements. 

The first, most illustrative, approach is the integral approach. This approach utilizes concrete models of 
thermoelements as thermoelectric power converters. For example, a model of a thermocouple element was 
employed by Ioffe. A model of thermomagnetic power converter was used by Harman and Honig. 
Such models are based on concrete variants of heat exchange between power converter, heat sources, and 
heatsinks, as well as the matching of thermoelement resistance with an external electric circuit wherein the 
mode of maximum efficiency or maximum electric power is attained. Under this approach the expressions 
for efficiency could be derived only for the few simplest cases. In general, the problem of efficiency 
determination is very complicated, and its general solution has not been obtained so far because of its 
awkwardness with account of temperature dependences of material properties, optimal configuration of 
thermoelements, and other peculiarities of thermoelements and their operating conditions. 

Therefore, the local approach is more productive than the integral one (Figure 2.1). 

This approach to efficiency calculation was proposed by Zener 16 and developed in Ref. [17] for 
thermomagnetic elements. Under this approach, parameters and properties of the thermoelectric energy 
conversion process in question are determined for the infinitesimal volume of thermoelement working 
medium. Then, if necessary, a change-over to integral parameters is made. 

In Ref. [14] this approach was used to create a generalized local thermoelement model suitable for 
materials with arbitrary anisotropy in the presence of a magnetic field and deformations. Let us consider 
this model in more detail. 

Figure 2.1 shows an infinitesimal volume dV of thermoelectric material. It is located between two 
isothermal surfaces with temperatures Tand T + dT. The lateral surface of volume dVis adiabatically 
isolated. A quantity of heat dQ h comes through the surface T + dT into the middle of the volume. Heat 
flow density on the surface T + dT is equal to q h . In the volume d there is thermal into electric energy 
conversion, resulting in electric current generation of density j. A quantity of heat dQ comes through the 
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surface T of volume dV! Thus, the spent quantity of heat is dQ h , and the resulting electric energy in 
conformity with the law of conservation of energy is dQ h — dQ c . Let us introduce the efficiency of 
thermoelectric energy conversion in volume d (local efficiency): 

i?iocai = dQ h - dQ c /dQ h (2.41) 

For its calculation it is advisable to introduce dimensionless parameters: dimensionless current i = 77* 
j//<*f, dimensionless electric field e = II"E/p" K*t, and dimensionless heat flow k = q//<*f, where t, j are 
modules of vectors VTand j. 

Then Ohm's and Fourier laws can be represented in a dimensionless form: 



e = Ri + Z TAto 
k = m - ©to 



(2.42) 
(2.43) 



with introduced tensors of material constants — dimensionless thermoemf A = a/a*, dimensionless 
resistivity]? = pip* , dimensionless thermal conductivity & = kIk* , and dimensionless Peltier coefficients 
tt= 77/77*. 

Expression (2.26) can be transformed as 



1?local = TJCarnotCP*/ + ajt/lftj ~ f 2 ) = 7J Carnot 5 \l + (UZ*T)H - 1) 



(2.44) 



where T/camot is the Carnot efficiency of a local model, and parameter 8 is determined by formula (2.11). 
Formula (2.29) takes into consideration a new parameter: 



Z*T= a*II*/K*p* 



(2.45) 



As can be seen from Equation 2.44, the local efficiency is a function of three factors, namely the Carnot 
efficiency, the operating mode of thermoelement determined by a dimensionless current i, and the value 
Z*T that for given temperature T is determined by the properties of thermoelement material and the 
magnetic field. 

Dependence Tji oca i on dimensionless current (Figure 2.3) has maximum t/™^ at 



i = i opt = Vi + z*r - i 



and the value of this maximum efficiency is 



max rj * 

^local — T ?Camot ze 



(2.46) 



(2.47) 
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FIGURE 2.3 Local efficiency 7ji oca] and coefficient of performance e loca] vs. dimensionless current for various Z*T 

values. 
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introducing a new parameter 



Ze* 



8(yf 



i/Vi 



i + z*t - i/vi + z*r+ 1 



(2.48) 



Let us elucidate physical content of parameters Ze* and Z*T. 

As is evident from Equation 2.47, parameter Ze* indicates the portion of the Carnot efficiency that can 
be obtained in local efficiency with the use of given thermoelectric material. Formula (2.36) implies that 
the generalized thermoelectric efficiency Z*T is the only parameter to determine Ze" \ hence the local 
efficiency for given material at temperature T. Ze as a function of Z" T is shown in Figure 2.4. 

From the figure it can be seen that at Z*T —> oo the efficiency of thermoelectric converter tends to the 
efficiency of the Carnot cycle. 

To go from the local efficiency T7i OC al to the integral efficiency T?i nt egrai> it is necessary in a thermoelectric 
energy converter to integrate the local efficiencies for all infinitesimal dV" volumes along current lines. 
This integration can be replaced by integration with respect to temperature T. 

Then 



max 
^/integral 



1 - exp(-S z ) 



(2.49) 



where 



i: 



Ze*dlnr 



(2.50) 



Graphically, the integral is depicted as area S z of a hatched curvilinear trapezoid under the curve of Ze" 
versus In T plot (Figure 2.5). The larger this area is, the higher is the integral efficiency of conversion. 
Parameters similar to Ze* and Z*T for partial cases of several thermoelement types are well known in 
the theory of thermoelements. Thus, in classical work by A.F. Ioffe, " from the consideration of an 
integral model of thermocouple thermoelement parameters were introduced to estimate the quality of 
materials used in thermocouple thermoelements, namely, the thermoelectric figure-of-merit of isotropic 
materials 



a uIk 



(2.51) 



having dimensions of the inverse temperature and a dimensionless parameter of thermoelectric figure- 
of-merit (Ioffe criterion) for such materials: 



Io = ZT 
The isotropic materials are characterized by quality criterion also introduced by Zener: 



(2.52) 



Ze = Vi+zr - l/Vl + ZT + 1 



(2.53) 
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FIGURE 2.4 Ze as a function of Z* T. 
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FIGURE 2.5 Calculation of integral efficiency as S 2 = J T 2 Ze d In T. 



It is readily apparent from Equations 2.33, 2.36, 2.37, and 2.38 that the well-known Ioffe and Zener 
criteria are partial cases of parameters Z"T and Ze" . Thus, for example, for the isotropic medium from 
Equation 2.32 we get the efficiency of thermocouple element 



17=17Carnot(M-l)/(M+T c /r h ) 



(2.54) 



where 



M = Vl + l/2(Z (T e + T h )) 



(2.55) 



Z = («i - a 2 f/(jK l p l + jK 2 p 2 ) 



(2.56) 



is an optimal thermoelectric figure-of-merit of thermoelement, a, p, and k with the respective indexes 
are thermoelectric coefficient, electric resistivity, and thermal conductivity, accordingly, of the first and 
second thermocouple legs. 

For the efficiency of the well-known anisotropic thermoelement from Equation 2.47 we also 
get the expression which coincides with Equation 2.54, but in this case M = y/l + Z ± T*, where 
Z ± — (Aa) 2 /4p* k ', Aa is the difference in components of the thermo-EMF tensor, and T* is the 
thermoelement temperature averaged along its gradient. 

For the Nernst thermoelement we have a similar expression, but the anisotropy of thermoemf is caused 
by a magnetic field and characterized by the value Q L B: 



Z ± = (Q ± B) z /p*K* 



(2.57) 



where Q 1 is the Nernst-Ettingshausen transverse coefficient. 

It is evident that parameters Z* T and Ze" comprise the well-known Ioffe and Zener criteria as partial 
cases. However, they are more general and suitable for an arbitrary anisotropy of thermoelement material 
and an arbitrary magnetic field. Therefore, it is advisable to call Z*T a generalized thermoelectric figure- 
of-merit, and accordingly, Ze" a generalized Zener criterion. 

The thermodynamic restrictions for the generalized criteria of efficiency and the Peltier and Seebeck 
coefficients are given in Table 2.2. 

The generalized expressions (Equation 2.47 and Equation 2.49) for the local and integral efficiencies 
are the most universal and allow calculating the efficiency of thermoelements for any anisotropic and 
inhomogeneous medium in an arbitrary magnetic field. 
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TABLE 2.2 Thermodynamic Restrictions for Thermoelectric Parameters 



Parameter 


Name 






Restriction 


Remark 


lo = ZT 


Ioffe criterion 






< ZT < oo 


For isotropic medium 


Ze 


Zener criterion 






< Ze < 1 




Ze' 


Generalized Zener criterion 






< Ze* < 1 


For media of any 


Z'"T 


Generalized thermoelectric figure 


■of- 


merit 


-1 <Z"T< oo 


symmetry and in the 


IT, a* 


Generalized Peltier and Seebeck coefficients 


IT > a*T 


presence of a magnetic 












field 



2.4 Conclusion 



Modern thermodynamic theory of thermoelectricity based on the thermodynamics of irreversible 
processes: 

(i) Gives new forms as given in Equation 2.31 and Equation 2.32 of the Thomson relations. 

(ii) Establishes new generalized quality criteria (Equation 2.45 and Equation 2.48) for material used 

in thermoelectricity, 
(hi) Discovers new thermoelectric effects (Table 2.1). 
(iv) Allows, based on the local expressions for efficiency (Equation 2.44) and temperature 

distribution, calculating the integral efficiencies (Equation 2.49) for any thermoelements 

(Figure 2.5). 
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